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Since metallic surface states on (111) noble metals are free-electron like, their propagators can 
be evaluated analytically. Since they are well-screened, one can use simple tight-binding formalism 
to study their effects. The needed phase shifts can be extracted from experiment. Hence, one can 
now make quantitative predictions of these slowly-decaying, oscillatory indirect interactions. For 
the (isotropic!) pair interactions (which decay as the inverse square of adatom-adatom separation), 
remarkable agreement has been obtained with experiments by two groups. We have extended the 
formalism to consider the full indirect ("triple") interaction of 3 adsorbates, which is the sum of 
the 3 constituent pair interactions plus the non-pairwise "trio" contribution, which tends to decay 
with the 5/2 power of perimeter. Here, we concentrate on interactions due to ordered overlayers 
and to linear defects, relating the latter to the interactions of (nxl) ordered overlayers and both to 
the constituent pair and trio interactions. We compare with experimental studies of interactions of 
adatoms with adchains and of consequent ID motion of adatoms trapped between two such parallel 
chains. We discuss implications for step-step interactions (on vicinal surfaces), with attention to 
the modification of the surface state itself for small terrace widths. 

Puisque le modele d'electron fibres s'applique pour les etats metalliques de surface (111) des 
metaux nobles, les propagateurs peuvent etre evalues analytiquement. Puisqu'ils sont bien ecrantes, 
on peut simplement employer le formalisme des liaisons fortes pour etudier leurs effets. Les 
dephasages necessaires peuvent etre extraits de l'experience. Par consequent, on peut maintenant 
faire des previsions quantitatives pour ces interactions indirectes oscillantes et decroissant lente- 
ment. Pour les interactions isotropiques de paires (qui decroissent comme l'inverse du carre de 
la distance entre adatomes), un accord remarquable a ete obtenu avec des experiences par deux 
groupes. Nous avons prolonge le formalisme pour considerer l'interaction indirecte de 3 adsorbants, 
qui est la somme des 3 interactions constitutives de paires plus la contribution a trois corps (trio), 
qui tend a decroitre comme la puissance de 5/2 du perimetre. Ici, nous nous concentrons sur les 
interactions dues aux couches adsorbees ordonnes et les defauts lineaires, reliant ce dernier aux 
interactions (nxl) des adorbats ordonnes et tous les deux aux interactions constitutives de paires 
et de trio. Nous comparons avec des etudes experimentales des interactions des adatomes avec des 
ad-chaines et du mouvement a ID des adatomes emprisonnes entre deux telles chaines parallele. 
Nous discutons des implications pour les interactions les marches (sur des surfaces vicinales), en 
faisant attention a la modification de l'etat de surface pour de petites largeurs de terrasse. 

PACS numbers: 73.20.Hb, 73.20.At, 68.35.Dv, 68.37.Ef 



INTRODUCTION AND PARAMETERS 

Metallic surface states, i.e. surface states crossing the 
Fermi level, have dramatic consequences that can be ex- 
plored at the atomic scale by modern surface probes such 
as scanning-tunneling microscopy (STM). Here, we sum- 
marize our progress to date in understanding the con- 
sequences of these states for nanoscale interactions not 
only between adsorbed atoms but also between chains of 
atoms and other atoms or chains. This work is prepara- 
tory to extensions to step interactions. We also present 
some helpful tabulations not published previously. 

While the simple asymptotic expressions for indirect 
interactions are valid only for separations larger than sev- 
eral/many atomic spacings, the more general expressions 
are valid for any surface-lattice separation (and could in- 
deed be generalized to arbitrary separations by allowing 
different phase factors for the interacting species). On 



the other hand, for atoms at nearest-neighbor spacings, 
particularly for homoepitaxy (or for adatoms larger than 
substrate atoms), direct interactions should also come 
into play, and are then expected to overwhelm any in- 
direct effects. Thus, for example, one must be cautious 
about using the formalism below to predict interactions 
between dimers or chains and atoms, specific any indirect 
interactions that involve a propagator between the two 
members of the dimer (or neighboring atoms in a chain). 

Table I summarizes key parameters that characterize 
the relevant isotropic Shockley surface states found on 
the (111) facet of the noble metals Cu and Ag. For 
both surfaces there now exist experimental investigations 
of the long-ranged adsorbate interactions 0, The 
table compares experimental values of the surface-state 
band parameters, obtained via STM with our large, 
sized-converged calculations by standard first-principles 
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TABLE I: Shockley surface-state parameters and Thomas- 
Fermi (bulk-screening) wavevectors of the Cu and Ag (111) 
surfaces. The Shockley band is characterized by the effec- 
tive electron mass m e ff, a Fermi energy £f (measured rel- 
ative to the bottom of the surface-state band), and a corre- 
sponding in-surface Fermi wavevector qF = ft -1 \ / 2m c geF and 
half wavelength Af/2 = tt/qf- The Thomas- Fermi screening 
lengths, kzh, are obtained as in Ref. Table adapted from 
refs. 0[Mlli0 





STM Cu(lll) 


DFT 


STM Ag(lll) 


DFT 


e F (eV) 


0.38 a 


0.42 


0.065 a 


0.045 


m of r/m e 


0.44 a 


0.38 


0.40 a 




q F (A" 1 ) 


0.21 a 


0.20 


0.083 a 




Af/2 (A) 


15.0 a 


15.5 


37.9° 




k T p (A) 


0.552 




0.588 




S adsorbate 


S F = ±tt/2 6 








Cu adsorbate 


S f = ±vr/2 c ' d 








Co adsorbate 


Sp = ±ir/2 d 




S F = n/3 d 





"Ref. [3] 'Ref. [6] 



R( :. [1] 



'Ref. [2] 



DFT the agreement is good. The table also shows 

STM measurements of the scattering phase shifts 5 F ^ 
reported 0,0 E| for various adsorbates, e.g. from stand- 
ing waves in "quantum corrals" 0. Finally the table 
shows estimates for the Thomas-Fermi screening length 
k FF 0- The surface-state electron response arises at 
much longer length scales, Af/2 3> fc^ F , and so will dom- 
inate the long-range adsorbate interaction. 



PAIR INTERACTIONS 

The interaction between adsorbates on a metal surface 
can involve an elastic, an electrostatic, and an indirect 
coupling through electronic states of the substrate. The 
long history of theoretical investigation of indirect adsor- 
bate interactions dates back nearly four decades 0; the 
history of this oscillatory, long range interaction has been 
amply documented ^cj- Lau and Kohn ^l| pointed out 
that the range of the interaction increases dramatically 
when the mediation is by a surface rather than a bulk 
states. Recent theory work 0, applied these ideas 
to the above-mentioned iso trop ic surface-state bands to 
find the pair- interaction fllll2|: 

A£ pair (d; S F ) = -Im / In (l - [t (e; S F )g Q (qd)} 2 ) (1) 
77 Jo 

a pasym,, , v ( 2sin(g F ) \ 2 sin(2 gF rf + 26 F ) 

~ A£ W W6 F ) = -6F ~ J—J^ -(2) 



( qF dy 



The simple analytic expression holds at asymptotic sep- 
aration d > Af/2. The effective T-matrix to(e;5 F ) — 
— (2h/m c s) sin(<5o(e)) exp(z<5o(e)), is determined by the 
s-wave phase shift <5o(e) with the boundary condition 



<5o(£f) = 5p. The surface propagator go{x) becomes ba- 
sically the cylindrical Hankel function of the first kind 



9o[x) 



meff (1) 

2h [x) 



,m e ff exp(ix — in/ 4) 



CO. 



(3) 

To obtain the simple asymptotic expressions, tg must be 
small enough so that ln[l — . . .] can be expanded to lead- 
ing order and x must be large enough to replace Hq 1 \x) 
by an outgoing circular wave. 

Subsequent STM measurements of Cu and Co adsor- 
bate dynamics on Cu(lll) and Ag(lll) 0,0 have ver- 
ified that the interaction has period Af/2 = n/q F and 
the quadratic decay of the envelope with separation, both 
without adjustable parameters. Accounting for the over- 
all magnitude requires insight into inelastic losses to bulk 
states. 



TRIO INTERACTIONS 

Study of the interaction of three adsorbates 0, Q| 
serves as a bridge from pair interactions to multi- 
adsorbate interactions in clusters. The three adsorbates 
are taken to bond to substrate positions i — 1,2,3. 
The triple-adsorbate cluster adsorption energy is calcu- 
lated [13| by combining a formal expansion p^fllla o f the 
adsorbate-cluster energy with scattering theory |12fl : 

A£ ; triplo(dl2, d,23, d 3 \\5 F ) 
3 

= AE p^T(dij;5 F ) + AE trio (di 2i d 23 ,d 3 i;S F ) (4) 
i>j=i 

= - Im / rfe In [l - {K 2 12 + K 2 3 + K 2 31 ) -2K 12 K 23 K 31 ] , 
n Jo 

where Kij is shorthand for to(e; S F )go(qdij). This 
triple-cluster interaction includes a new trio contribu- 
tion Ai? t rio which arises from constructive interference 
of electrons which traverse the entire cluster parame- 
ter G?i23 = di2 + d 23 + d 3 i. In the asymptotic limit, 
d-123 > 3Af, we obtain the analytical result |13lll4l |: 



4 f e " 

AE uio (d 12 ,d 23 ,d 31 ;S F ) ~ Im / deK 12 K 23 K 31 (5) 

k Jo 



-e F sin 3 (S F ) 



16V2~\ sm(q F d 123 + 3<5f - §7r) 
) 7123 (q F d 123 f/ 2 ' 



For completely absorbing scatterers the trio interaction 
result © is reduced by a factor of 1/8 (see Ref. [l3|). 
Since the scattering is taken to be s-wave, the trio in- 
teraction depends overwhelmingly on the perimeter d\ 2 z 
and is insensitive to the shape: the geometrical prefac- 
tor 7123 = y ^123/^12^23^31 varies little except for highly 
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TABLE II: Comparison of indirect interactions on surfaces 
mediated by [metallic] surface and bulk states and as well as 
bulk interactions (mediated by bulk states). The pair and 
trio decays refer to the envelope of the oscillatory interaction 
in the asymptotic regime. 





Surf, via surface 


Surf, via bulk 


Bulk 


A F /2 


~ 15.0A [Cu(lll)] 


~ 2.3A [Cu] 


~ 2.3A [Cu] 


Dispersion 


Isotropic 


Anisotropic 


Anisotropic 




e « (fifc||) 2 /2m* 




£n(k) 


Compu- 


Simple: para- 


Messy: multi- 


Messy 


tation 


bolic 2D band 


ple 3D bands 




Pair decay 


oc <r 2 


oc d~ 5 


oc d~ 3 




=£- observable 


=>■ insignificant 


RKKY 


Trio decay 


oc d~ 5 ' 2 


oc d~ 7 


ocd" 4 



distorted arrangements [l3l ll4T ] . Also > trio interac- 
tions can affect the barriers of atoms approaching; grow- 
ing clusters, an issue of recent theoretical study [lj| . 

Our results are summarized in Table 2. We emphasize 
that our calculations are non-perturbative, resulting in 
the physically-important phase shift Sp absent in pertur- 
bative approaches (e.g. Ref. 0). Since S F can differ for 
various adatom-substrate combinations (cf. Table 1), one 
can in principle select a system that will have a minimum 
at an arbitrary lattice spacing. 

While some evidence exists that the pair interaction 
alone is inadequate at non-asymptotic separations, there 
has not yet been a comparable experimental confirma- 
tion; trio interactions between adatoms and dimers are 
likely to be dwarfed by direct-interaction effects in the 
dimer, but other effects can be envisioned. 

The preceding process can be extended to compute in- 
teractions between 4, 5, and more adatoms. The formal- 
ism for bulk impurities, readily convertible to surfaces, 
was worked out by Harrison |17|. Alternatively, one can 
consider the interaction energy of superlattices of adsor- 
bates 0, 0] . In this way, one can relate the integrand 
for AE of a fractional overlayer to that of a full mono- 
layer. However, the simple expression for a full monolayer 
given in Ref. involves "tricks" related to the simple 
model employed that are subtle to generalize. 



INTERACTIONS WITH CHAINS 

By viewing a chain as the sum of its constituent atoms, 
one can readily add up these interactions to show 



A£ chain-atomW « " £ F Sin (S F ) — 



{q?£fl 2 



(6) 

where £ is the distance from the atom to the chain. The 
remarkable 3/2 power law was recognized over a decade 



ago [20j . A similar result should arise from consideration 
of the interaction energy of an (nxl) array of adatoms 
(but cf. warning at the end of section 3). 

Inserting parameters for Cu(lll) into Eqn. JJJJ, we find 
minima when I is 9, 24, 39, and 54 A. In counting the 
occurrences of atoms between 20 and 30 A from a chain, 
Repp |2l| did indeed find the behavior of Eqn. © . The 
chain- chain interaction has the same form as Eqn. J^) 
since the second chain can (also) be viewed as the sum 
of individual atoms, each of which have this interaction. 

An atom between two parallel chains will experience 
a ID corrugation potential parallel to the chains. Repp 
constructed such a situation for Cu atoms on Cu(lll) 
with atomic manipulation |2l| and produced STM movies 
of atoms wandering along the trough. Since the chains 
are of finite length, the well depth decreases near the 
ends of the chain. Hence, the atom is trapped in this 
furrow. If the chains are far enough apart, there are mul- 
tiple furrows. Repp 21] observed two atoms, in furrows 
nanometers apart, moving back and forth individually. 

One can imagine extensions of these ideas such as pro- 
ducing gridworks of chains with a regular set of traps 
for atoms or a maze of walls through which atoms might 
move as stupid rats. Computing the corresponding po- 
tential surface is then a fairly well-defined task. 



COMPLICATIONS IN GOING FROM CHAINS 
TO STEPS 

Surface states are not so robust as bulk states, so one 
cannot blithely view them as unaffected by the adsorp- 
tion process. Baumberger et al. [24| show that, on vicinal 
Cu(lll), the surface state is shifted up (and so q F re- 
duced) as the terrace width £ decreases. However, when 
the steps are decorated with CO, the energy shift be- 
comes downward with decreasing I! 

Furthermore, Ortega et al. [25j find that when £ de- 
creases sufficiently (in particular, when the misorienta- 
tion of Cu(lll) increases beyond 7°), the surface state 
is no longer that of the (111) facet but is determined by 
the vicinal surface itself. The periodic potential of the 
steps then opens a gap in the parabolic band structure. 

Both these arguments assume implicitly that the steps 
are straight and uniformly spaced, neither of which are 
generally true. It is not clear how the meandering of 
steps or the fluctuations in £ alter these results or, con- 
versely, how the interactions affect the meandering and 
distribution of the steps. (In concise words, are the steps 
"actors or spectators?" [26].) 

The existence of slowly-decaying oscillatory interac- 
tions should have profound implications for the distri- 
bution of terrace widths P(£). In general the dominant 
interaction between steps comes from entropic and elastic 
repulsions, both of which vary as £~ 2 . As a consequence 
P(£) has a "universal" form depending only on the ra- 
tio £/(t) (and the strength of the £~ 2 repulsion) but not 
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on the mean spacing (I), i.e., not on the misorientation. 
With surface states, this scaling breaks down, as has been 
observed experimentally |27| . 

Furthermore, the oscillatory interaction introduces a 
new length scale Xf- Thus, the equilibrium crystal shape, 
which is expected to be independent of crystal size, would 
seem to acquire some size-dependent behavior, at least 
for small crystallites. Since the £~ 3 / 2 decay of the enve- 
lope is slower than the £~ 2 of the pure repulsion, it is not 
clear what changes arise in the Pokrovsky-Talapov [2^] 
"critical behavior" of the curved regions near the edges 
of facets. 

CONCLUSIONS 

In summary, we present both an asymptotic evaluation 
and an exact model calculation for adsorbate interaction 
energies mediated by an isotropic Shockley surface-state 
band, as found on noble-metal (111) surfaces. While 
this interaction is primarily the sum of pair interactions, 
there can be significant trio corrections. Such interac- 
tions can pla y a role in the low-temperature adsorbate 
assembly ^3, 0, H3| , and efforts are being made to in- 
vestigate them directly 0, We can on this basis 
evaluate the interaction between a chain of adatoms and 
another chain and/or other adatoms. Novel nanostruc- 
tures can be imagined and actually contructed [2l| by 
skilled experimentalists. 

As noted, the slowly-decaying oscillatory interactions 
should affect a broad range of phenomena and should ap- 
ply to any situations in which defects create localized per- 
turbations on surfaces with surface states, e.g., magnetic 
interactions. Thus, the exchange coupling should oscil- 
late with the same period Af/2 as the adatom-adatom 
interaction; however, there is no a priori reason to expect 
that the phase shift Sf will be the same. Thus, one can 
imagine a rich phase diagram. 

In subsequent papers we will present a detailed investi- 
gation 0| of the surface-state-derived interactions asso- 
ciated with chains and nanostructures. We will also pro- 
duce a careful and thorough analysis and assessment of 
the assumptions involved in our approach 29] , with com- 
ments about extensions to systems in which, for example, 
rapid screening of the adsorption bond is questionable. 
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